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ABSTRACT: In this paper we introduce an approach to increase integration rate of elements of a dynamic 


comparator. In the framework of the approach we consider a heterostructure, which consists of a substrate 
and an epitaxial layer with special configuration. After that we consider doping by diffusion or ion 
implantation of several specific areas of the epitaxial layer with annealing of them in the framework of 
optimized scheme. The algorithm of manufacturing of the dynamic comparator leads to decreasing of 
dimensions of elements of the dynamic comparator with possibility to increase their density. We also 
introduce an analytical approach for prognosis of technological processes with account spatial and at the 
same time temporal variations of parameters as well as nonlinearity of the considered processes. 
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INTRODUCTION 

An actual and intensively solving problems of 
solid state electronics is increasing of integration 
rate of elements of integrated circuits (p-n- 
junctions, their systems et al) [1-8]. Increasing of 
the integration rate leads to necessity to decrease 
their dimensions. To decrease the dimensions are 
using several approaches. They are widely using 
laser and microwave types of annealing of infused 
dopants. These types of annealing are also widely 
used for annealing of radiation defects, generated 
during ion implantation [9-17]. Using the 
approaches gives a_ possibility to increase 
integration rate of elements of integrated circuits 
through inhomogeneity of technological 
parameters due to generating inhomogenous 
distribution of temperature. In this situation one 
can obtain by decreasing dimensions of elements 
of integrated circuits [18] with account Arrhenius 
law [1,3]. Another approach to manufacture 
elements of integrated circuits with smaller 
dimensions is doping of heterostructure by 
diffusion or ion implantation [1-3]. However in 
this case optimization of dopant and/or radiation 
defects is required [18]. 


In this paper, we consider a heterostructure. 
The heterostructure consist of a substrate and 
several epitaxial layers. Some sections have been 
manufactured in the epitaxial layers. Further we 
consider doping of these sections by diffusion or 
ion implantation. The doping gives a possibility to 
manufacture field-effect transistors framework a 


©2022 Mahendrapublications.com, All rights reserved 


dynamic comparator so as it is shown on Figs. 1. 
The manufacturing gives a possibility to increase 
density of elements of the dynamic comparator 
[4]. After the considered doping dopant and/or 
radiation defects should be annealed. Framework 
the paper we analyzed dynamics of redistribution 
of dopant and/or radiation defects during their 
annealing. We introduce an approach to decrease 
dimensions of the element. However it is 
necessary to complicate technological process. 


Fig. 1. The considered dynamic comparator [4] 
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METHOD OF SOLUTION 


In this section we determine spatio-temporal distributions of concentrations of infused and implanted 
dopants. To determine these distributions we calculate appropriate solutions of the second Fick's law 


[1,3,18] 
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Boundary and initial conditions for the equations are 
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The function C(x,y,z,t) describes the spatio- (including annealing), with changing 


temporal distribution of concentration of 
dopant; T is the temperature of annealing; Dc 
is the dopant diffusion coefficient. Value of 
dopant diffusion coefficient could be changed 
with changing materials of heterostructure, 
with changing temperature of materials 


D,.=D,(x, ya} I & e 


Here the function D, (x,y,z,T) describes the 
spatial (in heterostructure) and temperature 
(due to Arrhenius law) dependences of 
diffusion coefficient of dopant. The function P 
(x,y,z,T) describes the limit of solubility of 
dopant. Parameter vy e[1,3] describes average 
quantity of charged defects interacted with 
atom of dopant [20]. The function V (x,,z,t) 
describes the spatio-temporal distribution of 
concentration of 
Parameter V* describes the equilibrium 


radiation vacancies. 
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concentrations of dopant and_ radiation 
defects. We approximate dependences of 
dopant diffusion coefficient on parameters by 
the following relation with account results in 
Refs. [20-22]. 


* 2 


tere 
av; - GB) 
v) 
distribution of concentration of vacancies. 
The considered concentrational dependence 
of dopant diffusion coefficient has been 
described in details in [20]. It should be 
noted, that using diffusion type of doping did 
not generation radiation defects. In this 
situation Gi= G2= 0. We determine spatio- 
temporal distributions of concentrations of 
radiation defects by solving the following 
system of equations [21,22] 
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Boundary and initial conditions for these equations are 
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Here po =1,V. The function I (x,y,z,t) describes the 
spatio-temporal distribution of concentration of 
radiation interstitials; D,(x,y,z,T) are the diffusion 
coefficients of point radiation defects; terms 
V2(x,y,z,t) and /2(x,y,z,t) correspond to generation 
divacancies and diinterstitials; kiv(x,y,z,T) is the 
parameter of recombination of point radiation 


O® (x, y,z,t) O 


defects; kii(x,y,z,T) and kvv(xy,z,T) are the 
parameters of generation of simplest complexes of 
point radiation defects. 


Further we determine distributions in space and 
time of concentrations of divacancies @®/(x,y,z,t) 
and diinterstitials @(x,y,z,t) by solving the 
following system of equations [21,22] 
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Boundary and initial conditions for these equations are 
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D (x,y,Z,0)=far (xy,Z), Pv (x,y,Z,0)=fov (x,y,z). (7) 


Here Da,(x,y,z,T) are the diffusion coefficients of 
the above complexes of radiation defects; 
ki(x,y,z,T) and kv (x,y,z,T) are the parameters of 
decay of these complexes. 


We calculate distributions of concentrations of 
point radiation defects in space and time by 
recently elaborated approach [18]. The approach 
based on transformation of approximations of 
diffusion coefficients in the following form: 
D,(X,y,Z,T)=Dop[1+& gp{X,y,Z,T)], where Dop are the 
average values of diffusion coefficients, 0<¢,<1, 


lg-(x, y,z,T)|<1, p =LV. We also used analogous 
transformation of approximations of parameters 
of recombination of point defects and parameters 
of generation of their complexes: 
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Gv.v(X%y,Z,T)], where kopi,2 are the their average 
values, 0<é,v <1, O<e1 <1, O<evv<1, | giv(Xy,z,T)|<1, 
| gi(Xy,Z,T)|<1, \gv.v(xy,z,T)|<1. Let us introduce 
the following dimensionless variables: 
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We determine solutions of Eqs.(8) with conditions (9) framework recently introduced approach [18], i.e. as 


the power series 


i= 


P(x.7,6.9) = Ye, Le LOD, (7.77.6.9). (10) 


Substitution of the series (10) into Eqs.(8) and 
conditions (9) gives us possibility to obtain 
equations for initial-order approximations of 


concentration of point defects Loy %i19, 3) and 
Vio 27: ?, 9) and corrections for them 
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2>1. The equations are presented in the Appendix. 
Solutions of the equations could be obtained by 


standard Fourier approach [24,25]. The solutions 
are presented in the Appendix. 


Now we calculate distributions of concentrations 
of simplest complexes of point radiation defects in 
space and time. To determine the distributions we 
transform approximations of diffusion coefficients 
in the following form: Dao(x,y,z,T)=Doo,[1+ 
EG o/(X,Y,Z,T)], where Dogp are the average values 
of diffusion coefficients. In this situation the 
Eqs.(6) could be written as 
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Farther we determine solutions of above equations as the following power series 


Now we used the series (11) into Eqs.(6) and 


®,(x,¥,2.t)=Le,®,,(2y,z). (11) 
approaches [24,25] and presented in the 
Appendix. 


appropriate boundary and initial conditions. The 
using gives the possibility to obtain equations for 
initial-order approximations of concentrations of 
complexes of defects ®,0(x,y,z,t), corrections for 
them @,i(x,y,z,t) (for them i =1) and boundary and 
initial conditions for them. We remove equations 
and conditions to the Appendix. Solutions of the 
equations have been calculated by standard 


Now we calculate distribution of concentration 
of dopant in space and time by using the approach, 
which was used for analysis of radiation defects. 
To use the approach we consider following 
transformation of approximation of dopant 
diffusion coefficient: Di(xy,Z,T)=Do1[1+ 
é1g1(X,y,Z,T)], where Do. is the average value of 
dopant diffusion coefficient, O<er< 1, 
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Eq.(1) as the following series: 


C(uyat)= Le Le'C(ayzt): 


Using the relation into Eq.(1) and conditions (2) 
leads to obtaining equations for the functions 
Ci(x,y,z,t) (i 21, j 21), boundary and _ initial 
conditions for them. The equations are presented 
in the Appendix. Solutions of the equations have 
been calculated by standard approaches (see, for 
example, [24,25]). The solutions are presented in 
the Appendix. 


DISCUSSION 


In this section we analyzed spatio-temporal 
distributions of concentrations of dopants. Figs. 2 
shows typical spatial distributions _ of 
concentrations of dopants in neighborhood of 
interfaces of heterostructures. We calculate these 
distributions of concentrations of dopants under 
the following condition: value of dopant diffusion 
coefficient in doped area is larger, than value of 


1.0 


C(x.0) 


0.0 
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We analyzed distributions of concentrations of 
dopant and radiation defects in space and time 
analytically by using the second-order 
approximations on all parameters, which have 
been used in appropriate series. Usually the 
second-order approximations are enough good 
approximations to make qualitative analysis and 
to obtain quantitative results. All analytical results 
have been checked by numerical simulation. 


dopant diffusion coefficient in nearest areas. In 
this situation one can find increasing of 
compactness of field-effect transistors with 
increasing of homogeneity of distribution of 
concentration of dopant at one time. Changing 
relation between values of dopant diffusion 
coefficients leads to opposite result (see Figs. 3). 


Substrate 


Fig. 2a. Dependences of concentration of dopant, infused in heterostructure from Figs. 1, on coordinate in 
direction, which is perpendicular to interface between epitaxial layer substrate. Difference between values of 
dopant diffusion coefficient in layers of heterostructure increases with increasing of number of curves. Value of 
dopant diffusion coefficient in the epitaxial layer is larger, than value of dopant diffusion coefficient in the 
substrate 
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Fig. 2b. Dependences of concentration of dopant, implanted in heterostructure from Figs. 1, on coordinate in 
direction, which is perpendicular to interface between epitaxial layer substrate. Difference between values of 
dopant diffusion coefficient in layers of heterostructure increases with increasing of number of curves. Value of 
dopant diffusion coefficient in the epitaxial layer is larger, than value of dopant diffusion coefficient in the 
substrate. Curve 1 corresponds to homogenous sample and annealing time © = 0.0048 (Lx2+L,2+L,) /Do. Curve 2 
corresponds to homogenous sample and annealing time © = 0.0057 (Lx2+L,2+L,2)/Do. Curves 3 and 4 
correspond to heterostructure from Figs. 1; annealing times © = 0.0048 (L,2+L,2+L,?)/Do and © = 0.0057 (L:2 
+Ly2+Lz7) /Do, respectively 


It should be noted, that framework the considered 
approach one shall optimize annealing of dopant 
and/or radiation defects. To do the optimization 
we used recently introduced criterion [26-34]. The 


distribution by step-wise function yw (x,y,z) (see 
Figs. 4). Farther the required values of optimal 
annealing time have been calculated by 
minimization the following mean-squared error 


optimization based on approximation real 
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Fig.3a. Distributions of concentration of dopant, infused in average section of epitaxial layer of heterostructure 
from Figs. 1 in direction parallel to interface between epitaxial layer and substrate of heterostructure. 
Difference between values of dopant diffusion coefficients increases with increasing of number of curves. Value 
of dopant diffusion coefficient in this section is smaller, than value of dopant diffusion coefficient in nearest 
sections 
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Fig.3b. Calculated distributions of implanted dopant in epitaxial layers of heterostructure. Solid lines are 
spatial distributions of implanted dopant in system of two epitaxial layers. Dashed lines are spatial 
distributions of implanted dopant in one epitaxial layer. 
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We show optimal values of annealing time as 
functions of parameters on Figs. 5. It is known, 
that standard step of manufactured ion-doped 
structures is annealing of radiation defects. In the 
ideal case after finishing the annealing dopant 
achieves interface between layers of 
heterostructure. If the dopant has no enough time 
to achieve the interface, it is practicably to anneal 
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the dopant additionally. The Fig. 5b shows the 
described dependences of optimal values of 
additional annealing time for the same parameters 
as for Fig. 5a. Necessity to anneal radiation defects 
leads to smaller values of optimal annealing of 
implanted dopant in comparison with optimal 
annealing time of infused dopant. 


L 


x 


Fig.4a. Distributions of concentration of infused dopant in depth of heterostructure from Fig. 1 for different 
values of annealing time (curves 2-4) and idealized step-wise approximation (curve 1). Increasing of number of 


curve corresponds to increasing of annealing time 


C(x,0) 


Fig.4b. Distributions of concentration of implanted dopant in depth of heterostructure from Fig. 1 for different 
values of annealing time (curves 2-4) and idealized step-wise approximation (curve 1). Increasing of number of 


curve corresponds to increasing of annealing time 
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Fig. 5a. Dimensionless optimal annealing time of infused dopant as a function of several parameters. Curve 1 
describes the dependence of the annealing time on the relation a/L and = y= 0 for equal to each other values 
of dopant diffusion coefficient in all parts of heterostructure. Curve 2 describes the dependence of the 
annealing time on value of parameter ¢ for a/L=1/2 and € = y= 0. Curve 3 describes the dependence of the 
annealing time on value of parameter € for a/L=1/2 and ¢= y= 0. Curve 4 describes the dependence of the 
annealing time on value of parameter yfor a/L=1/2 and ¢= €=0 
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Fig.5b. Dimensionless optimal annealing time of implanted dopant as a function of several parameters. Curve 1 
describes the dependence of the annealing time on the relation a/L and €= y= 0 for equal to each other values 
of dopant diffusion coefficient in all parts of heterostructure. Curve 2 describes the dependence of the 
annealing time on value of parameter ¢ for a/L=1/2 and € = y= 0. Curve 3 describes the dependence of the 
annealing time on value of parameter € for a/L=1/2 and ¢= y= 0. Curve 4 describes the dependence of the 


annealing time on value of parameter yfor a/L=1/2 and e= €=0 


CONCLUSION 


In this paper we analyzed possibility to increase 
integration rate of elements of a dynamic 
comparator. To make the analysis we consider a 
heterostructure, which consists of a substrate and 
an epitaxial layer with special configuration. With 
account structure of the considered 
heterostructure we consider doping by diffusion 
or ion implantation of several specific areas of the 
epitaxial layer with annealing of them in the 
framework of optimized scheme. We obtained, 
that the algorithm of manufacturing of the 
dynamic comparator leads to decreasing of 
dimensions of elements of the dynamic 
comparator with possibility to increase their 
density. We also introduce an analytical approach 
for prognosis of technological processes with 
account spatial and at the same time temporal 
variations of parameters as well as nonlinearity of 
the considered processes. 
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APPENDIX 


Equations for the functions Ly (yon, —, 9) and Vie (wn, —, 9), i>0,j 20, k>0 and conditions for them 
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Solutions of the above equations could be written as 
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Equations for functions ®,i(x,y,z,t), i 20 to describe concentrations of simplest complexes of radiation 
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Equations for the functions Cj(x,y,z,t) (i 20, j =>0), boundary and initial conditions could be written as 


AC,,(x, y,Z,t) -D, a’C,,(x, y.Z,t) +D,, OC AR Zt) +D,, a’C,,(x, y, Zt) 
Ot ox Oy Ce 
AC,,(x,y,z,t) -p, OC. (x.y, ot) P a’C,,(x,y,z,t) Z OC, (x, y,z,t) 7 

Ot Ox Oy Oz 


Ox Ox 


PT ails v2.7) Coal 22) it 


Oy y 


F Dy, ae yer Sale d2) ,1>1; 


ce PC Ax YZt) &C,,(x, y,z,t) 
OL 2 OL 2 OL 2 
Ot Ox Oy Oz 


+D,, O C’ (x, y,z,t) OC,(x,y,z,t) pe. O C’ (x, y,z,t) OC, (x, y, zt) " 
Ax| P’(x,y,z,T) Ox Oy| P'(x,y,z,T) Oy 


OC\(xy, Zt) _ p, PCy(oy,z.t) 


+ 


Pankratov 


International Journal of Advanced Science and Engineering www.mahendrapublications.com 


Int. J. Adv. Sci. Eng. Vol.8 No.3 2288-2309 (2022) 2306 E-ISSN: 2349 5359; P-ISSN: 2454-9967 


8 | Cylx y, Zt) OCyg(x Zt) | 
* Az| P’(x,y,z,T) Oz , 


aC,,(x, y,Z,t) a. a? C,,(x, y, z,t) +Dy Oo Cle y,Z,t) +D,. a Calne y, z,t) 
Ot Ox Oy Ox 


+ 


Cr (x,y, zt) 
(x, y> it) P’(x, yz) 


x 


+ Dy, = C(x y, 2,8) (x, ¥,2t) OC 24) + e C,, 
Ox P(x, y,z,T) Ox dy 


4 ICH 9 Zst) + o C (x, y,z,t) 
dy Oz| 


Cz (xyz) 0Cylny.20)] 
P’ (x, y,z,T) Oz 


eee ree a 6 Eee 
OL 


P’(x,y,z,T) oz Ax| P’(x,y,z,7) 


x 


. C,,(x, 2 fe 0  ( y,z,t) 
dy Oz 


Culm yZt)|, 8 | Coleyzt) OC (yzt)|, 0 | Chlsyzt) OCalx yz) 
Ox dy| P’(x,y,z,T) Oy Oz| P’(x,y,z,T) Oz 


6C,,(x, y,z,t) D, °C, (, 924) , py PC, (%,y,2t) |p OC, (x.y, zt) 
i 2 L 2 


Ot Ox Oy ms Oz 


POR C. (x, y,2,t) Sia Zt) OC 9,20) le. 
Ox P’(x,y,z,T) Age ay 


x 


Caen) 
(x,y, it) P’(x, y,z,T) 


Lal Zt) |, 2 | (x,y,2,1) Sin, He drZet) OC ul YZ) a 
Oy Oz) °°” P*(x, y,z,T) Oz a 


O C’ (x, y,z,t) AC,,(x,y,z,t) O C’ (x, y,z,t) AC,,(x,y,z,t) 
oD + + 
Ax| P’(x,y,z,T) Ox Oy| P’(x,y,z,T) Oy 


_ 9 Genet) Cale .2) 1D, | oO er 27) 2Cul zz) e 


Og: P\sy.27) Oz Ox Ox 


rece g, (x, y,2,7)OLale 2-24) as O g, (x,y,z, OL r24) 
Oy Oy Oz az 


OC, (x,y,z) 2) AC, (x, y,z,t) -~0 BC, (x y,zt) 0 
Ox 4 Ox =e , Oy y=0 | 
Pankratov 


International Journal of Advanced Science and Engineering www.mahendrapublications.com 


Int. J. Adv. Sci. Eng. Vol.8 No.3 2288-2309 (2022) 2307 E-ISSN: 2349 5359; P-ISSN: 2454-9967 


AC(uyzt] 

Oy a | 
OC, (x,y,z) =0 A Cle ¥.t) = 0, 120, j>0; 
Oz aS v8 ee —_ 


Coo(x,y,Z,0)=fc (x,y,Z), Ci(x,y,z,0)=0, i 21, j 21. 


Functions Cj(x,y,z,t) (i 20, j >0) could be approximated by the following series during solutions of the above 
equations 


2 © 
Zi > F,.c,(xJe,(y)e,(zee(t). 


Col. ¥, 2st )= PS "TLL a 
eo ty 2 go ty 2 


Here Ec (t ) = exp 


Ly Ly 


E.= je,(u)fe,()f felu.va) c(w)dwdvdu; 


0 0 


20 : Ts 2 


Culeryoeat)=— ap eM Feclul Jey) ey Z)e ef ercl-F)L 8, (4) | c,(v)I g, (w.v,w.T) x 


n=l 0 0 0 


« ,(u) OD) gpa ydude-— > $nFe,(a)e,(y)e,(cJeu(Mfe(-2)* 


0 


. fe,(u) js, (rfc, (Je, (u,v, ws 7) Gaal wt), wdvdudt— Se Yn Felt) . 


XC, (x)e,(y)c, (z)fe,.(- r)fc. (u) fc, Ole (v)g, (u,v, wT) Cian wit) 4 wdvdudt,i 


0 0 0 
>1; 
t L Ly 


2 20 % : i 
Cale rot) =- ZF nkcelale.Q)oeldfed-afs. lhe W)fesbn) 


 CvalbeVWit)OC (UWE) ar a ae 2 
P’(u,v,w,T) Ou LL 


x y 


oe oe CU (uv wt) OCg(tv,¥, Qn 
x Jexcl- lesa) etree aly ue ") 4 wdvdudt— Ee yn e (t)x 


Fp et Ficu (wey y)ey()e elt) x 


<Fec.eJelMecDfesehal foyer Cal) graduate 


Pankratov 


International Journal of Advanced Science and Engineering www.mahendrapublications.com 


Int. J. Adv. Sci. Eng. Vol.8 No.3 2288-2309 (2022) 2308 E-ISSN: 2349 5359; P-ISSN: 2454-9967 


Wa t Ly Ly L, 
Cols yse8)=— eS nese vole Mel nnn 
x Cie oo (wv,wsz) OC ltevewst) Ween, rm F .c,(x)e,(y)x 


P’(u,v,w,T) Ou 


; om mn L; C’"(u,v,w,t) OC, (u,v, w, 
xn, (z)erclt)] enel- r)fe,(u)}s,(v)] C,,(u,v,w,T) ston altv w,T) . 


xc (w)dwdvdudt-— a 
LL 


TE Fcc a)e,(yey(2)e c(h euc(F) feu) fer) x 


xy Zz 


C*'(u,v,w,t) OC,,(u,v,w,7) 
P’(u,v,w,T) Ow 


L, 2m 
x s,(w)C,,(uv,w.) dwdvdudt— ELL, yn c,(x)x 


Fc, 9) e(2)eel*Megel—£)f5,(u) fe, (v)fc,(ov)C, (u,v, , 7) Sem) 


, : : Ou 
ye E00 (u, V,W,; Ja wdvdudr- 2 a nF oC. (xe, (y Je,(ze,-(t)fe,e(- r)[c,(u)x 
“Pru, V,W, T) ° ° 
xfs (v)fe (w)C aebidilelles a ¥ nx 
eee - P’ (u,v,w,T ) Ov LLL a 
Be A ey Sle C’"(u,v,w,T) 
x Fogle y)e, Verl Mterel MN. Cu) fest) fs, (w) Cy usv. 2} * 
0 0 0 0 P (u,v,w,T ) 
x Calo) a 4 ydudt-—*_$ F,c,(a)e,(y)e,(2)e lO) feuelt) 5,0) 
Ow LLL, a : ° 
«nf Ven i) Oleg wa yaude 7 Se,(s)eelt) 


Feed oNecl-ANfee) fs) fen) er eo aya vaude> 


LEE Fetsldes ve e)erleel-2) fee) feof sso)» 


: C’ (u,v,w,t) OC, (u,v,w,7) 
P’(u,v,w,T) Ow 


dwdvdudt; 


Pankratov 


International Journal of Advanced Science and Engineering www.mahendrapublications.com 


Int. J. Adv. Sci. Eng. Vol.8 No.3 2288-2309 (2022) 2309 E-ISSN: 2349 5359; P-ISSN: 2454-9967 


Ly 


2 00 t L, y L, 
C,(uy.2t)= 2 Sn co seve )elfevc(2) fs, (u) fes()featw)» 


0 
OC,,(u,v,w, 2n @ 
« 2, (un, w,7)2Calenner) LEE EM eeded eel)» 
x ne Ae Gis: Gyfe: (w)g, (u,v, wT )OLalte, wt) wdvdudt— ae x 
0 cone oa Av Lie 


XEN Cell lofe, AC t)fc Auf s(o)fs,( Jerlurm7) ee Dawa vdude x 


0 


x Fee ise (y)eu(@) Fp SE Petalades(v)ele)eclO eue(#) fs) elo) 


i C oC 2 es 
xn [c,(w) eae ae us ) 4 wdvdudt-— ae an F .c,(x)c,(y)x 


t Ly Ly L, C’ VW, OC, Vw, 
«cel fecl-eNfese) fs.) ers els) aya vad r 


“TBE hele aia eNle,tu)fe,()fs,) Sete), 


P’(u,v,w,T) 


OC, (uv, w,7) (xo(ye{z)e,.(fe,(-e)fs,(u)x 
Ow : : 
x fe (v)fe (w)C,,(wv.w tin r) OC luv Wit) 4 eine Ynx 
0” oO” eee P’(u,v,w,T ) Ou LE, ia 
x Fc,x)e,(y)6,(2)eeOMf eel) f,(u) f5,(v) fo, (or) ex ete e2) 2 Cole v.02) 
0 0 0 0 - (u,v, w,T) Ov 
x C, (u,v, w,t)d wdvdudt-— : vn F..c,(x)c,(y)c, (z)e,.(1) e,o(— r)fe, (u) ‘a 
7 n=l 0 0 


x je,(r)fs,(w)C (u,v, w,T) Cees Cale. Wt) 4 wdvdudt. 


All © 2022 are reserved by International Journal of Advanced Science and Engineering. This Journal is licensed 
under a Creative Commons Attribution-Non Commercial-ShareAlike 3.0 Unported License. 


Pankratov 


International Journal of Advanced Science and Engineering 


www.mahendrapublications.com 


